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Introduction  Research  and  testing  performed  with  the  equipment  purchased  under  this 
instrumentation  award  during  the  first  year  of  operation  consists  of  three  projects.  The  first 
was  an  identification  of  various  damping  mechanisms  in  a  damped  beam  made  of  quasi 
isotropic  composite  material.  The  second  was  an  identification  of  nonlinearities  in  a  joint 
used  in  the  ABES  structures.  The  third  was  in  conjunction  with  accelerometer 
measurement  used  to  verify  a  combined  velocity  and  acceleration  observation  method. 
Each  project  is  briefly  described  in  the  following: 

Damping  Identification  Several  intense  experiments  using  the  laser  vibrometer  were 
performed  to  measure  the  response  of  a  quasi  isotropic  cantilevered  beam  with  a  removable 
tip  mass  excited  by  an  impulse  at  various  locations.  This  data  was  collected,  stored,  and 
sent  to  AFOSR  researchers  at  Brown  University  sing  BITNET.  The  data  was  analyzed  in- 
reference  1  using  a  spline  based  estimation  procedure,  starting  with  a  partial  differential 
equation  model  of  the  structure.  A  clear  advantage  over  modal  methods  based  on  a  finite 
dimensional  model  of  the  same  system  was  observed*,  this  is  described  in  reference  1, 
which  is  attached.  ’  i 

Nonlinearity  In  response  to  a  query  from  the  Air  Force  Weapons  Lab,  a  single  bay  of 
the  ABES  structure  was  tested  using  the  laser  sensing  system.  The  goal  here  was  to 
identify  the  nonlinearities  that  occurred  in  this  structure  to  the  presence  of  the  joints.  It  was 
found  that  if  the  joints  are  properly  torqued,  there  exists  a  reasonable  level  of  amplitude  and 
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a  reasonable  frequency  range  for  which  the  system  behaves  in  a  linear  fashion.  This  is 
detailed  in  reference  2,  which  is  also  attached. 

Combined  Observation  A  two  signal  approach  to  observation  was  proposed  and 
tested.  This  method  uses  an  accelerometer  and  the  laser  system  to  simultaneously  measure 
both  velocity  and  acceleration.  With  this  information  it  is  illustrated  that  the  sampling 
theorem  can  be  violated  and  correct  modal  data  can  still  be  identified.  This  is  discussed  in 
reference  3. 

Control  Experiments  Several  control  experiments  are  currently  under  way  using  the 
laser  vibrometer  system.  However,  no  results  are  currently  available.  The  experiments 
under  way  are  1)  vibration  control  of  a  slewing  beam  using  laser  feedback  and  2)  structural 
control  of  a  truss  using  velocity  feedback  and  proofmass  actuators. 
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An  Experimental 

Investigation  of  the  Nonlinearities 
in  the  Joints  of  the  ABES 
Structure 


Introduction  A  single  bay  of  the  ABES  test  structure  was  experimentally 
investigated  for  the  purpose  of  clarifying  nonlinearities  thought  to  be  present  in  the 
joints.  The  joints  of  the  ABES  structure  are  a  threaded  male/female  configuration 
with  a  common  nut  used  to  lock  the  structure’s  threaded  members  into  the  joints. 
Joints  of  this  type  are  typical  to  most  large  space  structures  and  are  generally 
thought  to  be  a  source  of  non  linearity,  making  linear  system  identification  of  such 
structures  impossible. 

In  most  cases  it  is  desirable  to  be  able  to  use  linear  mathematical  models  for 
both  identification  and  control  of  a  structure.  Hence  the  joints  of  the  ABES 
structure  are  investigated  experimentally  for  nonlinear  behavior  in  hopes  of  finding 
and  classifying  regions  of  linear  behavior. 

Experimental  Proceedure  In  order  to  isolate  the  effects  of  the  joint,  a  single 
joint  with  to  member  was  configured  for  testing.  This  arrangement  was  suspended 
as  indicated  in  figure  1. _ _  _ _ 
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Experimental  Configuration  for  Joint  Tests 
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An  impact  hammer  (Kistler  model  number  9724A500  with  charge  amplifier  model 
number  5004)  was  used  as  input  to  the  test  structure.  Two  accelerometers  (Kistler 
model  number  8628 A5  with  charge  amplifier  model  number  5004)  were  used  to 
measure  the  resulting  acceleration.  The  output  signals  where  calibrated  against  and 
monitored  by  noncontacting  velocity  measurements  taken  by  a  DANTEC  Laser 
Vibrometer  System  (consisting  of  a  Hughes  Laser  model  number  3225H/PC,  a 
DANTEC  tracker  model  number  9055  No  214  and  a  DANTEC  frequency  shifter  model 
number  9055  No  012).  In  all  reported  data,  the  laser  velocity  measurements  agreed 
with  the  accelerometer  measurements.  A  GENRAD  2515  sixteen  channel  analyzer 
was  used  to  provide  the  data  and  calculate  the  appropriate  transfer  functions. 

Several  combinations  of  input  output  transfer  function  measurements  were 
made.  The  effect  of  impulse  amplitude  and  bolt  torque  (the  joint  tightness)  on  the 
resulting  non  linearity  of  the  frequency  response  was  examined. 

All  experiments  were  performed  in  the  Mechanical  System  Laboratory  at 
SUNYAB,  in  Buffalo,  New  York. 

Nonlinear  Behavior  The  measurement  of  and  parameter  identification  of  structures 
which  exhibit  nonlinear  behavior  is  complicated  by  several  effects.  First,  a 
fundamental  assumption  in  the  input/output  measurement  of  linear  structures  is  that 
the  transfer  function  between  any  two  points  on  the  structure  is  the  same  regardless 
of  which  point  is  chosen  as  the  input  and  which  point  is  chosen  as  the  output. 
This  assumption  is  used  extensively  in  standard  modal  testing  for  linear  structures  to 
reduce  the  number  of  calculations  needed  in  constructing  multiple  degree  of 
freedom  models  of  large  structures. 

A  second  nonlinear  phenomenon  that  leads  to  severe  error  in  linear 
measurement  and  identification  is  that  of  secondary  resonance.  The  common 
approach  to  structural  identification  in  linear  systems  is  to  measure  the  frequency 
response  of  the  test  structure.  The  peaks  that  show  up  in  the  magnitude  versus 
frequency  plots  are  interpreted  as  resulting  from  the  resonance  phenomena.  Hence 
each  peak  is  assumed  to  occur  at  a  natural  frequency  of  the  structure  and  the 
number  of  peaks  then  indicates  the  number  of  degrees  of  freedom  one  should  use 
in  a  lumped  mass  multiple  degree  of  freedom  model  of  the  test  structure.  In  fact, 
a  single  degree  of  freedom  linear  oscillator  equation  is  curve  fit  to  each  peak.  In 
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this  way,  the  parameter  estimation  yields  a  modal  damping  ratio  and  an  undamped 
natural  frequency  for  each  peak  assumed  to  be  a  single  degree  of  freedom  model. 
In  contrast,  the  force  response  of  a  single  degree  of  freedom  nonlinear  system 
produces  an  amplitude  versus  frequency  plot  with  multiple  frequencies.  This 
inherently  nonlinear  characteristic  is  refered  to  as  secondary  resonance  -  either 
superharmonic  of  subharmonic  depending  on  the  frequency  range  of  interest.  These 
secondary  resonances  cause  linear  parameter  estimation  methods  to  predict  incorrect 
parameters  as  well  as  a  false  model  order. 

To  illustrate  the  existence  of  subharmonics  and  superharmonics  in  the  testing 
of  structures  consider  a  simple  single  degree  of  freedom  model  of  the  structure  with 
mild  nonlinearity  in  its  stiffness  term.  A  well  studied  equation  is  Duffings’  equation, 
given  by 

X(t)  +  2euX(t)  +  u02X(t)  +  €  ccX3(t)  =  F0  cos  nt  (1) 

Here  X(t)  represents  the  displacement,  X(t)  the  velocity  and  X(t)  the  acceleration. 
The  quantities  e,  n  and  a  are  cons tr ants,  p  indicating  the  damping,  a  the  strength 
of  the  nonlinear  stiffness  and  €  a  scaling  factor.  The  driving  frequency  is  n,  the 
magnitude  of  the  forcing  input  is  F0  and  w0  represents  the  linearized  systems 
undamped  natural  frequency. 

The  solution  of  equation  (1)  can  be  obtained,  using  perturbation  methods,  for 
a  range  of  driving  frequencies.  First,  consider  the  steady  state  response  of  the 
system  of  equation  (1)  to  input  driving  frequencies  near  three  times  the  linear 
undamped  natural  frequency.  This  is  represented  by  n  =  3u0  +  eo,  where  a  is  a 
small  parameter  used  to  define  the  frequency  range  of  interest. 

The  response  of  (1)  can  be  shown  to  be  of  the  form 

X(t)  =  a(t)  cosgt  -  J-XO]  -  §—5  cos  nt  ♦  0(0  (2) 

Here,  0(«)  represents  all  those  terms  of  order  e  and  higher,  and  the 
functions  a(t)  and  7(t)  vary  slowly  with  time.  The  frequency  content  of  X(t)  reveals 
oscillation  at  two  frequencies  rather  than  one  as  in  the  linear  case.  The 
appearance  of  the  second  frequency  of  oscillation  is  called  a  subharmonic  oscillation 


because  it  occrus  at  1/3  of  the  driving  frequency.  A  standard  modal  analysis  of 
such  a  system  would  imply  a  two  degree  of  freedom  system  with  natural  frequencies 
of  u0  and  Uq/3,  rather  than  a  nonlinear  single  degree  of  freedom  system. 

Likewise,  examination  of  the  solution  of  the  steady  state  response  of  (1) 
around  the  driving  frequency  3fl  =  w0  +  eo  yields  a  solution  of  the  form 

X(t)  =  a(t)  cos  [3nt+7(t)l  +  ^ ^  cos  nt  +  0(0  (3) 

ft  -Ho 

In  this  case  a  superharmonic  appears  in  the  solution.  An  identification  of  the 
frequency  content  of  this  signal  would  contain  the  natural  frequency  and  the 
superharmonic  frequency  and  again  yield  an  erroneous  linear  result. 

The  following  experimental  results  represent  a  first  attempt  to  clarify  the  linear 
range  with  both  frequency  and  amplitude  of  excitation  for  the  joints  in  the  ABES 
structure. 

Results  As  mentioned  above,  one  clear  sign  of  the  presence  of  nonlinearities 
is  the  violation  of  the  condition  for  linear  transfer  function  that  the  frequency 
content  of  the  transfer  function  with  point  A  as  input  and  point  B  is  output  should 
be  the  same  as  the  transfer  function  with  point  B  as  input  and  point  A  is  output. 
The  top  plot  of  Figure  2  illustrates  the  magnitude  of  the  transfer  function  with  A 
as  input  and  B  as  output.  The  top  plot  of  Figure  3  is  the  magnitude  of  the 
transfer  function  with  B  as  input  and  A  as  output.  The  frequency  content  in  the 
lower  frequency  range  of  0  to  900  HZ  appear  the  same.  However,  in  the  frequency 
range  above  1000  HZ  the  two  spectra  are  quite  different  indicating  that  the  joint  is 
nonlinear  (joint  torqued  to  2  ft.  lbs.)  in  the  higher  range  with  excitation  frequencies. 

The  next  set  of  tests  examine  the  effect  of  input  amplitude  on  the  joint  transfer 
function  near  a  primary  resonance  (about  255  HZ)  with  the  joint  torqued  to  2  ft.  lbs. 
The  amplitude  of  the  impulse  to  the  joint  increases  starting  with  figure  4  and  moving  to 
figure  6.  Note  that  at  some  (not  measured)  level  of  input  amplitude  super  and  sub 
harmonics  occur.  This  is  clearest  in  figure  6.  The  existence  of  super  and  sub 
harmonics  is  of  course  a  clear  indicator  of  nonlinear  behavior. 


Figure  7  illustrates  the  same  frequence  range  and  fundamental  harmonic  as 
figures  4,  5  and  6.  However  in  this  case,  the  bolts  in  the  joints  are  only  tightened 
to  1  ft.  lb.  of  toque  instead  of  2  ft.  lb.  Again,  this  configuration  produces  a 
variety  of  super  and  sub  harmonics  indicating  nonlinear  behavior. 

The  third  set  of  experiments  tests  the  joint  linearity  for  high  levels  of  joint 
torque  (35  ft.  lbs.).  Figure  8  illustrates  a  low  amplitude  hammer  hit  and  figure  9 
illustrates  a  high  amplitute  hammer  hit.  The  resonant  frequencies  for  the  high 
torque  joint  remain  the  same  indicating  the  for  high  torques  the  ABES  structure 
behaves  in  a  linear  fashion. 

Summary  The  joints  for  the  ABES  structure  are  definitely  capable  of 
illustrating  nonlinear  behavior.  The  nonlinearities  in  the  joint  are  a  function  of 
three  factors  1)  input  amplitude  2)  torque  or  joint  tightness  and  3)  frequency  range. 
If  it  is  desired  to  model  the  ABES  structure  with  a  linear  model  it  is  necessary  to 
use  evenly  torqued  joints  of  about  35  ft.  lbs.  and  to  realize  that  at  some  level 
*  (high  amplitude)  of  disturbance  excitation  the  linear  model  will  fail  for  high 
frequency  excitations. 

It  is  the  conclusion  of  this  study  that  the  ABES  structure  can  be  effectively 
modelled  with  a  linear  model  for  lower  amplitude  disturbance  provided  all  the 
joints  are  properly  torqued. 
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PARAMETER  IDENTIFICATION  TECHNIQUES  FOR  THE  ESTIMATION 
OF  DAMPING  IN  FLEXIBLE  STRUCTURE  EXPERIMENTS 


H.T.  Banks  and  Y.  Wang 
Center  for  Control  Sciences 
Division  of  Applied  Mathematics 
Brown  University 
Providence,  RI  02912 

Abstract 

We  report  on  our  use  of  spline  based  inverse 
procedures  to  estimate  damping  coefficients  in  distributed 
parameter  systems  for  flexible  structures.  Damping  models 
involving  viscous  (air)  damping  and  Kelvin-Voigt  damping 
in  an  Euter-Bernoulti  framework  are  used  to  analyze  data 
from  vibration  experiments  with  composite  material  beams. 

I.  Introduction 

In  this  presentation  we  report  on  continuing  efforts 
that  we  have  pursued  for  the  past  year.  These  efforts 
involve  a  combination  of  experimental  investigations  with 
the  simultaneous  development  of  the  mathematical  and 
computational  aspects  of  a  theoretical  methodology  to 
support  the  experiments. 

Our  long  term  goals  include  the  understanding  of 
damping  mechanisms  in  complex  distributed  structures 
constructed  from  composite  materials.  Our  quest  has 
obvious  motivation  from  and  potential  consequences  for  the 
design  of  control  (active  and  passive)  systems  for  large 
flexible  space  structures. 

The  initial  efforts  on  which  we  report  here  entail  the 
testing  aed  development  of  spline  based  approximation 
techniques  for  inverse  problems  arising  in  attempts  to 
model  and  quantify  damping  in  composite  material  beams. 
Substantial  previous  mathematical  efforts  along  with 
numerical  tests  on  simulated  data  can  be  found  in  the 
literature  (see  (BCR],  [BR]  for  a  sample  of  some  of  these 
results  and  references  to  other  work). 

We  briefly  outline  the  underlying  mathematical  ideas, 
describe  the  experimental  equipment  employed,  and  then 
report  on  two  representative  experiments  and  our  findings 
using  our  techniques  with  the  data  from  these  experiments. 


The  basic  mathematical  model  that  we  have 
considered  in  connection  with  the  efforts  on  which  we 
report  in  this  note  consists  of  the  equations  for  a 
cantilevered  Euler-Bernoulli  equation  with  tip  mass, 
damping,  and  an  applied  transverse  force.  Specifically  we 
take  as  our  model  the  system  for  transverse  vibrations  of  a 
long  slender  (length  1)  cantilevered  beam: 

'  &  * EI  &  fe  *  <*•■> 

-  f(t.xp.T),  0  <  x  <  1,  t  >  0, 

«(«.°)  -  -fjHt.O)  -  0.  t  >  0  (2.2) 

EI  to}  (t,l)  +  ‘o1  a^8T(M)  “  °’  ‘  ”  °*  (2.3) 

El  (t.l)  ♦  cDI  (t,l)  -  m  (t.l).  t  >  0,  (2.4) 


u(0.x)  -  *(x),  (0,x)  -  <«x),  0  <  x  <  I.  (2.5) 

Here  p  is  the  linear  mass  density,  m  is  the  tip  mass,  cDI 
and  7  are  damping  coefficients  for  Kelvin-Voigt  and 
viscous  (air)  damping,  respectively,  and  f  is  a  parameter 
dependent  (xp  for  the  point  of  impact,  t  for  the  duration) 
term  representing  a  transverse  force  such  as  that  obtained 
with  a  force  hammer.  The  functions  p  anJ  <t>  represent 
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given  initial  displacement  and  velocity  of  the  beam. 

We  assume  that  we  are  given  observations  -  in  this 
case,  of  tip  acceleration  -  from  which  we  wish  to  estimate 
unknown  parameters.  In  general  these  include  the  stiffness 
coefficient  EI  as  well  as  the  damping  parameters  cDI,  y 
and  the  input  parameter  T.  Denoting  this  collection  of 
unknown  parameters  by  q  -  (EI,cDI,7,r),  we  may  formulate 
this  inverse  problem  as  the  least-squares  estimation  problem 
of  minimizing  the  functional 

J(Q>  “  ?  I  ft*1'  (*,,,:q)  *  (2  6) 

over  an  admissible  class  Q  of  parameter  values.  In  this 
fit-to-data  criterion  (a(t())  are  the  observations  (data)  for 
the  tip  acceleration  (at  x  -  I)  of  the  beam. 

This  estimation  problem  cannot,  of  course,  be  solved 
analytically:  however,  one  can  employ  iterative  optimization 
schemes  coupled  with  an  approximation  method  for  the 
infinite  dimensional  syttem  (2.1)  -  (2.5).  For  the  efforts 
discussed  here,  we  approximated  the  system  via  Galerkin 
procedures  using  cubic  B-spline  elements  as  explained  in 
(BGRW].  This  results  in  an  approximate  solution  uN  which 
satisfies  a  readily-solved  finite  dimensional  system 
approximating  (2.1)  •  (2.5).  We  use  this  in  the  criterion 
(2.6)  in  place  of  u,  thereby  obtaining  a  family  of 
approximating  estimation  problems  involving  minimization 
over  Q  of 

jN(q)  “  ?  i&n (*>,i:q)  ’  ,(,i)l* 

Solving  these,  we  obtain  a  sequence  of  estimates  (q*  for 
best-fit  parameters  in  the  original  estimation  problem  A 
theory  for  convergence  of  these  estimates  can  be  obtained 
(see  e.g.  (BCRUBR)JBGRW]  for  details),  but  we  shall  i.ot 
discuss  this  aspect  of  the  computational  techniques  lee 
since  the  focus  of  this  report  is  the  US£  of  th  ■  t 
techniques  in  studying  damping  mechanisms  in  composi  < 
material  structures. 

We  have  assumed  here  that  the  observations  are  takei 
with  tip  accelerometers  since  that  is  precisely  the  case  for 
the  experiments  discussed  below.  However,  we  could 
equally  well  assume  velocity  observations  are  taken;  indeed, 
we  have  also  used  these  methods  with  experiments  in 
which  a  laser  vibrOmeter  was  used  to  obtain  velocity 
measurements.  These  and  related  findings  will  be  detailed 
in  a  separate  manuscript  currently  in  preparation. 

Furthermore,  in  our  findings  discussed  here,  we  have 
concentrated  on  attempts  to  characterize  the  damping  in 
terms  of  Kelvin-Voigt  and/or  viscous  (air)  damping.  There 
is  considerable  evidence  to  suggest  that  damping 
mechanisms  in  composite  materials  are  significantly  more 
complex  than  the  ones  described  by  a  combination  of 
viscous  and  Kelvin-Voigt  damping.  Indeed,  current 
conjectures  include  models  that  involve  hysteretic  or 
hereditary  (time  and/or  space)  effects.  Recent  theoretical 
results  (BIJ  provide  a  convenient  mathematical  framework 
for  treating  inverse  problems  for  such  models;  we  are 
currently  using  the  associated  computational  techniques  to 
investigate  the  possibility  of  describing  experimental  data 
in  terms  of  models  with  hysteretic  damping. 


To  test  the  above  described  estimation  procedures  and 
initiate  our  damping  mechanism  studies  a  series  of 
experiments  were  carried  out  at  the  Mechanical  Systems 
Laboratory  located  at  SUNYAB.  Cantilevered  composite 
beams  with  a  removable  tip  mass  were  used  as  the  test 
structures  Each  beam  consisted  of  a  one  meter  long  piece 
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of  pultruded  quasi-isotropic  composite  used  by  DOD  and 
NASA  in  studies  of  potential  materials  for  the  proposed 
Space  Station  (WMJ.  The  beams  rectangular  cross  sections 
measured  on  the  order  of  6  mm  by  100  mm  provide  a 
lower  cost,  higher  specific  strength  alternative  to  extruded 
aluminum  alloy  structures. 

t  _  The  experiments  consisted  of  recording  the  vibrations 
after  a  free  release  from  an  initial  bent  configuration  or 
after  an  impulse  was  imparted  to  a  resting  beam  using  a 
Kistler  Instrument  Corporation  force  hammer,  instrumented 
with  a  load  cell  to  provde  measurements  of  the  input  force 
magnitude.  Kistler  accelerometers  (Model  8628  piezo  beam) 
were  attached  to  the  free  tip  of  the  beams.  The  data 
acquisition  was  performed  with  a  GENRAD  2515  analyzer. 
Time  response  data  was  collected  and  sent  via  B1TNET  to 
Brown  University  for  use  in  the  above  described 
algorithms.  A  copy  of  data  was  also  used  as  input  to  a 
standard  modal  analysis  package  [SDRC]  residing  in  the 
GENRAD  as  well  as  to  an  Eigensystem  Realization 
Algorithm  residing  in  a  VAX78S. 

The  modal  analysis  results  were  then  interpreted  in 
terms  of  the  partial  differential  equation  model  given  by 
equations  (2.1)  -  (2.5)  above.  The  measured  natural 
frequencies  were  used  as  data  in  a  simple  average  to 
produce  initial  estimates  for  E.  The  damping  parameters 
were  calculated  by  using  an.  appropriate  number  of  modal 
damping  ratios.  This  procedure,  detailed  below,  represents 
a  relatively  standard  engineering  approach  to  estimating  the 
parameters  in  a  vibration  test  of  a  distributed  mass 
structure.  These  parameters  were  then  used  in  a  numerical 
simulation  of  the  response  and  compared  with  the  response 
predicted  by  the  spline  based  methods  presented  here.  Note, 
that  this  modal  analysis  procedure,  however,  requires  that 
closed  form  solutions  of  the  frequency  equations  be 
available.  This  is  a  shortcoming  not  shared  with  the 
spline  based  procedures. 

IV.  Results 

We  present  our  findings  and  some  preliminary 
conclusions  on  two  representative  examples.  Further 
discussions  and  results  from  other  experiments  will  be 
given  in  a  paper  currently  in  preparation. 

Example  FR1.TM86 

In  this  experiment  we  used  a  composite  material 
beam  of  length  Im,  width  ,102m,  thickness  ,006m  with 
collocated  tip  mass  and  accelerometer.  The  linear  mass 
density  was  p  -  .021546  slug/ft  and  the  tip  mass  was  m  - 
.057182  slug.  The  beam  was  cantilevered  at  x  «  0;  it  was 
displaced  from  equilibrium  (with  this  initial  displacement 
measured  for  use  as  initial  data).  At  time  t  <*  0,  the 
beam  was  released  and  the  free  oscillations  were  then 
observed;  2048  time  observations  at  At  -  .007813  seconds 
intervals  were  recorded.  Thus  time  domain  data  over  the 
period  (0,16]  seconds  were  available.  Upon  taking  the  FFT 
of  this  data,  we  found  that  only  one  mode  was  observable 
in  the  oscillations;  this  was  at  1.875  Hz.  For  our  model 
we  used  (2.1)  -  (£5)  with  El  calculated  with  experimental 
modal  information  to  be  126.7083  lb  ft1. 

In  our  initial  parameter  estimation  efforts,  we  tried 
to  fit  the  data  over  (0,2]  seconds  with  a  model  featuring 
only  viscous  damping  (i.e.  CqI  -  0).  With  our  spline 
inverse  procedures  we  found  7  -  .0315  slug/(ft  sec)  and 
obtained  the  model  response  depicted  in  Figure  1. 
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■/'  . 

We  also  used  modal  techniques  to  try  to  recover  a 
value  for  7  in  (2.1)  -  (2.5).  That  is,  we  used  standard 
modal  techniques  to  estimate  {,  the  %  critical  damping  on 
the  mode  excited,  as  l  »  .269.  Since  one  can  readily 

argue  that  7  -  2(/Elp/»,<  this  approach  yields  the 
calculated  value  7ctl  -  .0826.  Using  this  value  in  the 
model,  one  does  not  obtain  a  very  good  fit  to  the  data. 
This  is  shown  in  Figure  2. 
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Figure  2 


However,  these  initial  findings  are  somewhat 
misleading.  Even  the  spline  inverse  procedures  don't  yield 
a  good  fit  of  the  model  (2.1)  •  (2.5)  with  only  viscous 
damping  if  we  consider  the  time  domain  response  over  a 
long  period  of  time.  The  best  results  we  could  obtain  are 
depicted  in  Figure  3  where  we  note  we  also  optimize  over 
m  (this  was  done  to  adjust  the  m  value  since  the  "tip 
mass*  was  actually  a  small  tip  body  centered  at  x  -  .963m 
or  3.7cm  from  the  tip  end).  The  ’optimal*  value  of  7  was 
7  -  .015174. 
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Several  important  features  of  our  findings  here  arc 
worth  noting.  (i)  This  is  a  "1  mode  -  1  damping 
parameter*  model  fit  in  either  case.  Our  findings  reveal 
that  such  a  situation  does  not  automatically  lead  to  success. 
If  the  wrong  damping  model  is  hypothesized,  even  good 
quantiative  methods  won't  yield  satisfactory  results,  (ii)  It 
is  important  to  use  the  data  over  a  long  time  interval 
whether  we  are  using  FFT  data  or  time  domain  data. 

Example  1R3.1NTMI186 

Again  a  composite  material  beam  of  length  1  m  was 
used.  The  width  of  this  beam  was  .104m,  the  thickness 
.00524m.  No  tip  mass  was  attached  but  as  usual  a  tip 
accelerometer  was  available  for  data  collection.  We 
observed  the  response  of  the  beam  after  it  was  hit  with  an 
impulse  hammer  (soft  tip)  at  xp  -  ,7m.  Three  dominant 
modes  were  observed  in  the  response,  at  3.7Hz,  23.3Hz  and 
61  Hz,  respectively.  The  beam  parameters  were  p  «  .02144 
slug/ft  and  El  -  107.103  lb-fr  (calculated  approximately 
using  modal  information).  We  first  optimized  on  El  to 
obtain  a  frequency  match  between  our  model  and  the  FFT 

of  the  data.  A  value  of  El  -  109.2367  lb -ft*  was 
obtained.  Modeling  of  the  input  (f  in  equation  (2.1))  is,  in 
general,  nontrivial  since  it  is  not  strictly  an  impulse.  We 
used  the  input  signal  (from  a  transducer  in  the  hammer)  to 
model  the  shape,  height,  and  duration  of  the  input  to  the 
beam.  We  then '  used  our  spline  inve^je  procedures  to 
estimate  the  two  damping  parameters  7  and  c„I,  obtaining 

values  7  -  .004235  slug/(ft  sec)  and  cDI  -  .005385  slug 
fts/sec.  The  corresponding  time  domain  response  for  the 
model  and  the  FFT  are  shown  in  Figures  5  and  6.  As  is 
obvious,  we  obtained  a  very  good  fit  to  the  data  (with  3 
modes  excited)  using  only  two  damping  mechanisms  in  (2.1) 
-  (2.5). 


Figure  3 


If.  however,  we  hypothesize  that  the  damping  is  Kelvin- 
Voigt  only  (e.g.  7-0,  Cj,I  nonzero),  our  procedures  give  a 
much  better  fit  to  the  data.  We  obtained  cpl  -  .1338  with 
the  time  response  and  corresponding  FFTs  shown  in 
Figure  4. 
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Figure  4 
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Figure  6 


If  we  try  to  estimate  7  and  c„I  using  modal 
techniques  (via  %  critical  damping  and  principles  relating  i 
to  7  and  CqI  similar  to  that  noted  in  the  discussion  of  the 
previous  experiment),  we  obtain  values  7„,  -  .2982  slug/(ft 
sec)  and  Cjjl^  -  .451 1  slug  ft*/sec.  If  we  use  these  in 
the  model  equations  (2.1)  -  (2.5),  we  find  that  the  response 
is  damped  out  much  too  quickly  and  does  not  resemble 
the  experimental  data. 

A  final  comment  on  our  findings  in  this  example  has 
important  implications  for  use  of  modal  methods  in 
attempts  to  understand  damping  mechanisms  in  distributed 
parameter  structures.  While  we  found  that  the  viscous 
damping  (7)  was  extremely  important  in  the  overall 
modeling  of  the  low  frequency  mode  and  the  Kelvin-Voigt 
damping  (CpI)  was  most  important  in  modeling  correctly 
the  two  high  frequency  modes,  the  damping  does  not 
decouple  on  the  modes.  We  belive  this  poses  serious 
difficulties  for  any  attempts  to  understand  damping 
mechanisms  using  only  modal  techniques. 
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